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Canonical structure of a generalized time-periodic harmonic oscillator is studied by finding the 
exact action variable (invariant). Hannay's angle is defined if closed curves of constant action 
variables return to the same curves in phase space after a time evolution. The condition for the 
existence of Hannay's angle turns out to be identical to that for the existence of a complete set of 
(quasi)periodic wave functions. Hannay's angle is calculated, and it is shown that Berry's relation 
of semiclassical origin on geometric phase and Hannay's angle is exact for the cases considered. 
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Berry's phase the geometric part of a change in the 
phase of a wave function under a cyclic evolution, has at- 
tracted great interest both theoretically and experimen- 
tally . A significant generalization of the Berry's phase, 
relaxing the adiabatic approximation, has been given by 
Aharonov and Anandan Q . The price we have to pay for 
this generalization is that the quantum states do not nec- 
essarily return to the original states up to a phase after 
the evolution of Hamiltonian's cycle, so that the geomet- 
ric phase may not be defined through the method in . 
This geometric phase may have a natural classical corre- 
spondent: Hannay's angle Q . Hannay's angle or geomet- 
ric phase has long been studied for the (generalized) har- 
monic oscillator Q . In a recent study |^ of the oscillators 
with time-periodic parameters, it has been proven that, 
if two linearly independent homogeneous classical solu- 
tions of the oscillator are bounded for all the time, there 
always exists a complete set of (quasi)periodic wave func- 



tions for an oscillator without driving force. For driven 
case, the ratio of the period of periodic particular solu- 
tion to that of Hamiltonian should be rational for the 
existence of the complete set. This oscillator provides 
an ideal system for the study of dynamics of Gaussian 
wave packets which has been applied to many problems 
in atomic and molecular physics as well as in quantum 
optics 0. 

Recently, the Hannay's angle for generalized harmonic 
oscillators and its relation to geometric phase were stud- 
ied by Ge and Child §, and by Liu et al. [|. Though 
Hannay's angle was originally formulated relying on the 
existence of action variable in the canonical structure of 
a system ^,|| , and the exact invariant (action variable) 
for the harmonic oscillator of time-dependent frequency 
has been given by Lewis more than thirty years ago, 
none of the works for Hannay's angle of the oscillator has 
benefited from the exact action variable. 



In this Letter, we will study Hannay's angle of the os- 
cillator by generalizing the exact action variable, and will 
compare the angle with the geometric phase. The model 
we will consider is described by the Hamiltonian: 
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The smooth functions M{t), w{t) and F{t) denote the 
positive mass, real frequency and external force, respec- 
tively, while the overdots denote differentiations with re- 
spect to time. In the Lagrangian description, the terms of 
smooth functions a{t), b(t) and f{t) can be interpreted as 
the results of adding total derivative terms to Lagrangian 
which do not affect the classical equation of motion . 
We require a periodicity for every coefhcient, so that 



H{x,p, t + t) — H{x,p, t). 



(3) 



The classical equation of motion for the Hamiltonian is 
given as 



^(Mi) + M{t)w^{t)x = F{t). 



(4) 



The general solution of this equation is a linear combi- 
nation of a particular solution Xp{t) and two linearly in- 
dependent homogeneous solutions u{t),v{t). As in pT| , 
u{t) and v{t) are dimensionless, while Xp{t) has the di- 
mension of length. For later use, we define p{t) and a 
time-constant as 



p= ^u^(^t) +v^{t), n = M{t)[v{t)u{t) -u{t)v{t)]. (5) 



The p satisfies the following differential equation 



(6) 



For the case that M is constant and a,b,F, f are 
zero, the Hamiltonian describes the harmonic oscillator 
of time-dependent frequency w{t). In this case, the exact 
invariant 
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whose level surfaces form ellipses in the phase space of 
any time t, has been found in |]l0| (For a derivation of 
Iq from the Hamiltonian of a simple harmonic oscillator, 
see [|l2|). The subscript is to denote that variables are 
defined when M is the unit mass and a,b, F, f are zero. 
The area enclosed by the ellipse is given by 27r/o, as can 
be easily proven by making use of Stokes' theorem 
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with a parameterization 
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Though the ellipse evolves as time passes due 



to the time dependence of po? the area enclosed by the 
curve or /q is a time-constant. 

To find action variable for the Hamiltonian in (1), we 
make use of the unitary transformations given in pT[ | . By 
applying the unitary transformations to /q, one can find 
that the transformed operator / is written as 



1 il'^ 

+ {{Mp + 2Map)(x ~ xp) - p{p - pp)}% (9) 

where pp is defined as pp = Mxp + 2Maxp + b. The gen- 
eralized harmonic oscihator is a system where the path 
integral for the kernel is Gaussian. Making use of this 
fact, the quantum theory of the model has been studied 
in Ref. [0. In this paper, we will consider the I{x,p,t) 
in (|9|) as a classical object. 

One can explicitly check that / is a constant of motion 
satisfying 



dl _ dl 
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where the subscript PB denotes that the term is a Pois- 
son bracket. Eq. (10) shows that / is the action variable 
of the oscillator system described by the Hamiltonian in 
(1); The value of / for a point of phase space at a given 
time stays constant along the trajectory generated by H. 
Again, by making use of Stokes' theorem as in (8), one 
can show that 2-kI is the area of ellipse of constant / in 
the phase space. The term with coefficient a changes the 
shape of ellipse, while the terms with coefficients h, F 
just move the center of the ellipse and the purely time- 
dependent term in the Hamiltonian has no effect on the 
ellipse. We parameterize the ellipse as 



cosQ — 
sin Q = 
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[I,tanQ]pB = -— tanQ 
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shows that Q is the angle variable. For the evalua- 
tion of p as a function of / and x, we should con- 
sider the two branches of an ellipse divided by the line 
p — (M^ -|- 2Ma){x — Xp) + Pp. In the branch above 
(below) the line, the momentum is written as 



p = Pp + (M- + 2Ma){x - Xp) 



1 

±-J2m-—{x~xp)-^, 
p V p 



(13) 



with upper (lower) sign. A generating function of the 
canonical transformation from {x, p} to action-angle vari- 
ables is given in the branch above (below) the line as 



F2{x, I,t) ^ S + J f{z)dz + MaXp 4- hxp + MxpXp 
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with upper (lower) sign. In (p_4|), 5 is defined through the 
relation 
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One can find that p = gives the relation (O) and 



Q = is compatible with the parameterization (|l 
The Hamiltonian for the same system in terms of action- 
angle variables is given as 



H = H{xiQ,I,t),piQ,I,t),t) + —^ 
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showing again that / is the action variable. It is note- 
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worthy that is a single-valued function of Q as 
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Hannay's angle is defined when the closed curves of 
constant action variables return to the original curves af- 
ter a time evolution [Hj^. For the generalized harmonic 
oscillator, this can be satisfied if both of p{t) and Xp{t) 
are periodic with a period which is an integral multiple of 
T. In the cases where such period exists, we shall denote 
r' as the period. The evolution generated by H{x,p,t) 
transports a family of ellipses through the phase space. If 
r' exists, then the given family of such ellipses get trans- 
ported back to itself after a time t' . Hannay's angle is 
defined as the integral of angle-averaged value of Q — ^ 
for the period r': 

1 pto+r' p2Tr o 0771 

Q^-^l I ^(^U^.(Q..))« (18) 

where to is an arbitrary time. After some algebra, one 
can find, making use of (6), that the Hannay's angle is 
written as: 



1 r' 
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A general generating function ^2(0:, /, t) for the action 
variable I{x,p,t) of (^ may be written as F2{x,I,t) = 
F2{x, I, t) + IQc{t) + 5{t). Qc{t) must be dimensionless, 
while 5{t) has the physical dimension of 5{t). With this 



general generating function, the angle variable Q is given 
through the parameterization of (|l^) where Q is replaced 
by Q — Qc{t), so that 

21 



(20) 



—pcos{Q - Qc) + Xp, 

^r{— + 2Ma)pcos(Q - Q^) 
\l p 



sin((3 - Qc) +Pp- 



This parameterization is compatible with Q = ^§f-, and 
implies that Qc{t) represents the rotational motion of the 
line in phase space from which the angle variable is mea- 
sured. The Hamiltonian in terms of / and Q has two 
additional terms to that of (16). Making use of the def- 
inition dl^), the Hannay's angle Qh for the generating 
function F2 is evaluated as Qh = Qh + A(5(to), where 
AQ(io) = Qcih + t') - Qc{to). Formally, Qc{t) can be 
any function of t, and thus Qh depends on to unless Qc(t) 
is periodic with the period r'. The condition for the ex- 
istence of meaningful Hannay's angle, therefore, is that 
the Qc{t) is periodic with the period r'. The condition 
is satisfied if the mapping from {x,p) to {Q,I) is periodic 
with the period r'; If the mapping has the periodicity, 
the value of Q(x,p, tg + r) is equal to that of Q{x,p, to) 
in addition to the time-periodicity of I{x,p,t) discussed 
above. From now on we only consider the case of the pe- 
riodic Qc{t), so that Qh does not depend on Qc{t) and 
is equal to Qh- In (19), one can find that there could 
be cases where the p{t) and thus the Qh depend on the 
way of choosing u{t) and v{t), since the mapping is deter- 



mined by the choice of classical solutions; However, for a 
given mapping of the periodicity, the Hannay's angle is 
unique. 

In general cases, it looks like that there is no reason of 
Hannay's angle being described by an integral of canoni- 
cal variables (see Ref. ) . For the generalized harmonic 
oscillator considered here, the Hannay's angle happens 
to be written as 



The expression of Hannay's angle in ( |21| ) may be facil- 
itated in easily finding the fact that the angle does not 
depend on Xp and pp] The terms containing Xp and pp in 
the integral of ( pT] ) which comes from the linear terms in 
the Hamiltonian are removed through the differentiation 
with respect to / and the angle-average. The Hannay's 
angle does not depend on the motion of center of the 
ellipses, reflecting the fact that the angle is a geometric 
quantity. If there is no linear term in the Hamiltonian so 
that Xp — Pp — 0, averaging over a half range of the an- 
gle variable which corresponds to a branch of the ellipse 
divided by any straight line passing through origin gives 
the same Hannay's angle. 

In the quantum treatment of the Hamiltonian in (1), 
a set of wave functions satisfying Schrodinger equation is 
given as pH 



X exp\-\M(t)a{t)x'^ + (M{t)xJt) + b(t))x\ 
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For > 0, iprn is square-integrable. The set of wave 
functions {'iprn{x,t)\ m — 0, 1,2,---} are complete with 
a choice of two linearly independent solutions u{t),v{t), 
and a particular solution Xp{t) |[Tl|], while, with the same 
choice of 17 > 0, the mapping from {x,p) to {Q,I) at a 
given time is one-to-one. 

The condition for the existence of a complete set of 
(quasi)periodic wave functions is satisfied, if p{t) and 
Xp{t) are periodic with a period which is an integral mul- 
tiple of r. Therefore, the condition of the existence of 
a family of time-periodic closed curves in phase space 
needed for the definition of Hannay's angle is exactly 
same to that of the complete set of (quasi)periodicity. 
In 1^, it has been shown that, if a homogeneous solution 
diverges then such complete set does not exist and thus 
Hannay's angle can not be defined. Making use of Flo- 
quet's theorem, it has also been proven that [^, if there 
exist u(<), v{t) finite all over the time, which are the cases 
we will consider from now on, there always exist two lin- 
early independent homogeneous solutions which give pe- 
riodic p{t) with period r or 2t; For the case of Xp — 0, 
there exists a complete set of (quasi)periodicity and thus 
a family of closed curves of periodicity. For the case of 



Xp ^ 0, the ratio of the period of periodic Xp to r should 
be rational for the existence of Hannay's angle and the 
complete set(s) of (quasi)periodicity. 

If t', an integral multiple of t, is the common period of 
periodic p and Xp, the overall phase change of the wave 
function ipm under the r'-evolution is given as 

x^--i^+l)[ + (23) 

The expectation value of the Hamiltonian for a wave 
function ^pm{x,t) is given as 
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Geometric phase for the wave function Tp„iix,t) is thus 
written as 
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In the Hannay's angle, the effects of the linear terms 
in the Hamiltonian has been removed through the angle- 
average and differentiation with respect to /. In the cal- 
culation of geometric phase, the expectation values of the 
linear terms can not depend on the m due to the orthog- 
onality of Hermitian polynomials and their recurrence 
relations. A simple relation is thus satisfied between the 
geometric phase and the Hannay's angle 



In Ref. 1^, the relation (26) has been suggested for gen- 
eral models; However, the suggestion was made through 
a semiclassical treatment of geometric phase and using, 
in fact, the Eq.(|2l|) as a defining relation for Hannay's an- 
gle which may not be equal to the definition of Eq.(p^) 
in general cases. A similar relation for the case without 
linear terms is given in Ref. 

In summary, we analyze the canonical structure of a 
generalized harmonic oscillator by finding the exact ac- 
tion variable. The analyses has then been used to find 
Hannay's angle which is defined relying on the time- 
periodic closed curves of constant action variables in 
phase space. Hannay's angle for the model considered 
can be defined if and only if there exist a complete set 
of (quasi)periodic wave functions whose geometric phases 
have been given in this paper. There could be cases where 
both of the Hannay's angle and the geometric phase de- 
pend on the way of choosing classical homogeneous solu- 
tions, while the angle is unique for a given choice of clas- 
sical solutions of the periodicity. It should be of interest 
if similar analyses would be possible for other models. 
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